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Abstract 

We have calculated jet shapes in low Q 2 ep and pp collisions in perturbation theory 
at order for the hard parton-parton processes. For the 'yp process resolved and direct 
contributions are superimposed. The dependence of the jet shapes on transverse energy, 
rapidity, and inner cone extension is studied. The numerical results of the calculation 
are compared with recent data from ZEUS at HERA and from CDF and DO at the 
TEVATRON. Good agreement is achieved if the problem of merging overlapping jets is 
taken into account by varying the parameter R sep as a function of transverse energy and 
rapidity. 
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1 Introduction 


A large fraction of the final state in hadron-hadron, electron-hadron (both high Q 2 and low Q 2 , 
the latter also being referred to as photon-hadron), and electron-positron collisions consists of 
high energy jets. These jets have an extended structure which can be studied experimentally 
and theoretically. A possible measure of this structure is the jet shape or jet profile, which 
depends on variables like transverse energy and rapidity of the jets, the jet algorithm, and the 
extension of jets. The jet shape measure is the function p(r, R , Et, p) where Et is the transverse 
energy of the jet and p is its rapidity with </>, its azimuthal angle, integrated out. This function 
p measures the fractional Et profile, i.e. given a jet sample with transverse energy Et defined 
with a cone radius R , p{r, R, Et,p) is the average fraction of the jet’s transverse energy that 
lies inside an inner cone with radius r < R. 

This jet shape p for the special value r = 0.2 has been measured for the first time in pp 
collisions at a/s = 630 GeV by the UAl collaboration [1] at the SppS collider. More detailed 
measurements of p in pp collisions at a/s = 1800 GeV at the TEVATRON have been performed 
with the collider detector CDF |2|. Here, the dependences of the jet shape p on the transverse 
energy Et and on r were studied for jets produced in the central region p] and compared 
with theoretical calculations. The measured jets at CDF have sufficiently high Et of order 100 
GeV, so that one can assume that gluon emission effects are much more important than long¬ 
distance fragmentation processes or soft interactions with “spectator” partons in determining 
the jet shape function p(r, R, Et,p)- This means that at sufficiently high energies, the shape 
of the jet should be calculable by perturbative QCD alone ignoring fragmentation effects and 
other soft interactions. Such calculations of p based on af finite-order perturbative QCD for 
pp collisions were first performed by Ellis, Kunszt and Soper (EKS) 0 and the results were 
compared to the CDF data J2|. 

The production of high transverse energy jets by quasi-real photons on protons (low Q 2 
ep scattering) has some similarity to jet production in pp collisions. As is well known, two 
mechanisms contribute to the photoproduction of jets at large transverse energy. The photon 
can either interact directly or via its quark and gluon content, the resolved contribution, with 
the partons originating from the proton. The resolved photoproduction dominates at HERA 
for lower Et and positive rapidities p (we assign p > 0 in the direction of the incoming proton). 
Near the maximum of the rapidity distribution (p ~ 2), the resolved and direct parts become 
comparable at Et — 30 GeV. Due to the dominance of the resolved component, jet shapes in 
low Q 2 ep collisions should be rather similar to jet shapes in pp jet production. Next-to-leading 
order (NLO) QCD calculations of jet shapes at HERA for the resolved component have been 
presented some time ago ||. Their similarity to jet shapes in pp collisions has been studied 
in || showing that due to a simple scaling behavior jet shapes for different c.m. energies and 
reactions can be correlated. This scaling behavior was used by the HI collaboration to compare 
jet profiles as measured in pp reactions || with the jet profiles obtained from pp data [l]]. 

Recently, the ZEUS collaboration presented their first results for jet shape functions p in ep 
collisions 0. They studied jet shapes p(r, R, Et,p) in inclusive single-jet and dijet photopro¬ 
duction for jets with E T > 14 GeV, — 1 < p < 2, and R = 1 as a function of E T , p, and the inner 
cone radius r. These measurements show that p changes appreciably with p and Et for fixed 
r. In this work we shall calculate p for the experimental conditions of the ZEUS measurements 
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in NLO including direct and resolved contributions to find out whether the behavior of p with 
varying p and E? found by ZEUS can be explained by fixed-order perturbative calculations. In 
particular we shall investigate how the merging of partons in relation to the merging procedure 
in the experimental analysis can be exploited to obtain a satisfactory description of the exper¬ 
imental results. With the same attitude we shall investigate jet shapes in pp collisions at the 
TEVATRON and compare our results with measurements of CDF [f|] and DO || in different 
rapidity and transverse energy regions. 

The outline of the paper is as follows. In section 2 we discuss in detail the merging procedure 
in the theoretical calculations and in the experimental analysis of low Q 2 ep collisions. In section 
3 we compare our numerical results with the recent ZEUS data. The pp jet shapes are considered 
in section 4. The results of the calculation and the comparison with CDF and DO data are 
given. Section 5 summarizes our results and ends with some concluding remarks. 


2 Jet Shapes in Low Q 2 ep Collisions 

2.1 Theoretical Input 

We consider positrons with energy E e = 27.5 GeV colliding with protons of energy E p = 820 
GeV. The photon emission is described by the Weizsacker-Williams equivalent photon approx¬ 
imation with a formula as in i with Qmax = 4 GeV 2 and 0.2 < y < 0.85, where y = E 1 /E e 
is the fraction of the initial positron energy transferred to the photon, in accordance with the 
conditions of the ZEUS measurements 0. For the resolved contribution we need the parton 
densities of the photon. We take them from the work of Gluck, Reya, and Vogt m converted 
to the MS scheme. The parton distributions of the proton are taken from CTEQ4M in the 
MS scheme [11]. Both are NLO parametrizations. The details of these structure functions are 


not essential, since p is a ratio of cross sections. For a s (p) we employ the one-loop formula 
with Nf = 4 and with A taken from the proton density (A^ = 296 MeV). We decided on 
the one-loop expression for a s since the jet shape is calculated only to lowest non-trivial order 
0{a 8 ). The renormalization scale p is set equal to the factorization scales and is put equal to E T . 


For calculating the jet profile p we need a jet definition. We adopt the cone algorithm 
of the Snowmass convention 0. According to this definition a jet is defined as transverse 
energy E? deposited in a cone of radius R in the rapidity-azimuthal angle plane. The jet axis 
is determined from 

W = X] Et^/Et, (1) 


z G cone 


<t>j — ^2 ET,i4>i/E T 

i G cone 

with the transverse energy calculated from 

Et = E T ,i■ 

i G cone 

A parton with kinematic variables (■, </>,) is included in the jet if the condition 

(Vi ~ Vj ) 2 + (& - </> j ) 2 < R 


( 2 ) 


(3) 


(4) 
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is satisfied. Since we have only up to three partons in the final state, not more than two partons 
can be combined into one jet. Since the jet shape p is fully determined from the 2^3 parton 
transition cross section we have the situation that the final state consists either of three or 
of two jets, where one of the jets resulted from the recombination out of two partons in the 
3-parton final state. 


As is well known this jet definition has two problems: double counting and merging. The 
condition (4) is equivalent to 

ETi + Et) 


\j(vi ~ vj) 2 + 0* - <Pj) 2 < 


ma x(E Ti ,E T ) 


R, 


(5) 


i.e. if the parton angles of partons i and j satisfy this condition, they are counted as one 
combined jet with transverse energy Et + Etj, but not as two smaller jets with energies Et, 
and Et ■ For the case that 


R < ( Vi - Vj) 2 + ((pi - 0j) 2 < 


Et + Ej 


ma x(E Ti , E t 


;R, 


( 6 ) 


which occurs, for example, for E T , = Et = Et/2, the two partons i and j might with some 
justification count also as separate jets. Therefore one has the choice to count only the com¬ 
bined jet {i and j} or to count {*}, {j}, and {i and j} as separate jets which would contribute 
to the inclusive single-jet cross section. We follow EKS |l3j and count only the combined jet in 
this case. This way we avoid any double counting of jets. In practice it is found, however, that 
the jet shape p is unaffected by the double counting issue of jet recombination except at r ~ 0 


The issue of merging is much more severe. It is related to the problem whether the cone 
algorithm used in the theoretical NLO calculation matches the jet definition in the experimen¬ 
tal analysis. This problem cannot be solved easily since in the experimental measurement of 
p one starts with the observation of hadrons which are recombined into protojets with a cone 
algorithm similar to the Snowmass definition. These protojets are recombined further until 
the conditions of the algorithm which include merging of overlapping cones are satisfied. The 
details of the algorithm used by the ZEUS collaboration will be described later. It is clear, 
however, that without further investigations it is not possible to decide that cone algorithms 
used for three parton final states match the cone algorithms used in the recombination of final 
state hadrons in the ZEUS analysis. 


The theoretical jet algorithm we are using will merge two partons into a single jet whenever 
the condition (5) is satisfied. This includes also the configuration when two partons with equal 
transverse energy Et, = E = Et/2 are just 2 R apart. In the calculation this is counted as one 
jet of transverse energy Et with a cone centered between the two partons. This means that for 
this case the Et of the jet is distributed at the edge at the cone and not in the vicinity of the 
center. Whether such a configuration will be counted as a single jet also in the experimental 
analysis will, of course, depend on the details of the jet algorithm used in the experiment. This 
problem was recognized by EKS 0 when they calculated jet shapes p for pp collisions and 
compared them with the CDF data J2J. To simulate the merging and other conditions of the 
CDF cone algorithm in a simple way they added the extra constraint, that two partons, i and 
j, when separated by more than R sep < 2 R, i.e. when 

\!(Vi ~ Vj) 2 + 0* - 0j) 2 > Rsep, (7) 


3 








are no longer merged into a single jet. With this additional constraint the fraction of E r near 
the edge of the cone is reduced. EKS also found that R sep = 1.3 R describes the jet shape p 
as a function of r for Ex = 100 GeV jets as measured by CDF very well. This means that for 
simulating the CDF cone algorithm the R sep constraint with R sep = 1.3 R should be introduced 
on the parton level. Of course, this value for R sep accounts only for the ideal situation of high 
Ex jets at CDF with very small fragmentation, underlying event, and other effects, which might 
distort the relation of the CDF cone algorithm to the algorithm employed on the parton level. 
From this we conclude that the jet merging in the CDF analysis can be simulated on the parton 
level with the parameter R sep < 2 R which reduces the merging of partons at the cone edge. For 
E t — 100 GeV jets R sep ~ 1.3 R. For lower E T ’s we expect additional broadening of the jets 
which we might simulate with R sep values slightly larger than 1.3 R. Therefore we consider R sep 
a parameter which is not fixed completely and depends on the kinematical parameters of the 
final state jets. Before we decide which R sep values should be assumed for the photoproduction 
of jets at HERA we shall describe the jet algorithm employed in the ZEUS analysis. 


2.2 Jet Algorithm in the ZEUS Analysis 


In the ZEUS jet search and the measurement of the jet shape p [7fl a cone algorithm JT2[] is used 
to construct jets on the basis of the energy depositions in the uranium-scintillator calorimeter 
(CAL) cells in both data and simulated events and also from the final state hadrons in the 
simulated events. Whether hadrons or CAL cells are considered is not essential. The following 
steps are taken to construct jets. CAL cells are considered with their p and (j) determined from 
the unit vectors joining the vertex of the interaction and the geometric centers of the cells. In 
the first step, each CAL cell with a transverse energy above 0.3 GeV is considered as a seed for 
the search. These seeds are combined if their distance in p — 0 space R = \J (A p) 2 + (A0) 2 < 1. 
The cone radius R = 1 is drawn around each seed and the CAL cells within that cone are 
combined to form a cluster. The axis of the cluster is defined according to the Snowmass con¬ 
ventions [P~2f . This yields ?/ c i us ter and Cluster which follow from the transverse energy weighted 
mean pseudorapidity and azimuthal angle of all the CAL cells belonging to this cluster (see 
(1), (2), and (3)). A new cone of radius R = 1 is then drawn around the axis of the resulting 
cluster. All cells with their geometric center inside the cone are used to recalculate a new cluster 
axis. The procedure is iterated until the content of the cluster is stable. Up to this point the 
algorithm is very similar to the jet definition at the parton level which we use to produce our 
predictions for comparison with the experimental data. Of course, with only three partons in 
the final state, the procedure is much simpler. Iterations are not needed and the possibilities 
for merging are very limited, since maximally only two partons can be merged. 


The difference to the theoretical algorithm occurs in the treatment of overlapping clusters. 
In the ZEUS analysis two clusters are merged when the common transverse energy exceeds 
75% of the total transverse energy of the cluster with the smallest transverse energy. Otherwise 
two different clusters are formed and the overlapping cells are assigned to the nearest cluster. 
Finally a cluster is called a jet if the corrected Ex exceeds 14 GeV. Concerning the rapidity 
only jets with p in the range — 1 < p < 2 are selected. 

We emphasize that in this experimental procedure of defining jets there are several steps that 
have no analogy in the theoretical jet definition with only three partons in the final state. First, 
in the experimental definition iterations are needed to define stable clusters. Second, there is 
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the merging of two clusters depending on the shared energy. In the theoretical definition there 
is neither an iteration nor shared energy. The latter would only occur in a theoretical analysis 
with at least four partons in the final state, i.e. if we went to one order higher in the QCD 
calculation. In this case we have NNLO two-jet, NLO three-jet, or LO four-jet production. 
In O(o.g), where there is only NLO two-jet and LO three-jet production, merging of partons 
occurs only once, so that iterations and overlapping cones do not arise. 


The energy sharing and merging of overlapping clusters in the ZEUS analysis is the same 
procedure as used by the CDF collaboration to define jets in pp collisions. In both measure¬ 
ments, jets are constructed with the choice R — 1. Therefore we expect that the jet definition in 
the ZEUS analysis corresponds to R sep ~ 1.3i? in the merging of partons for the theoretical jet 
algorithm at least for the ideal case of large Ej- and/or jet production in the backward direction. 
For negative rf s the jet production is a superposition of nearly equal contributions from direct 
and resolved production with x 7 ~ 1 (x 7 is the fraction of the initial positron energy trans¬ 
ferred to the photon). In this region there is little disturbance of the jets through additional 
interaction with the remnants. In the region x 1 > 0.75 the measured transverse energy flow 
around the jet axis is described reasonably well without extra multiple interactions included in 
TT| . This is already the case for rather low Et but improves the larger the 

parameter into the theoretical jet definition can 


the Monte Carlo 
Et is. Of course, the introduction of the R sep 


only be a phenomenological device to model the much more involved experimental jet definition. 


Further information on R sep can be gained from a study of experimental and theoretical 
jet definitions for photoproduction |T3[. In this study inclusive dijet cross sections da/dp with 
E t > 6 GeV and \rji — r] 2 \ < 0.5 for direct and resolved photoproduction as a function of the 
average pseudorapidity of two observed jets rj = (rji +r] 2 )/2 were compared for three jet defini¬ 
tions denoted EUCELL, PUCELL and KTCLUS. The first two definitions are cone algorithms 
whereas KTCLUS is a cluster algorithm. PUCELL is identical to the jet definition described 
above. KTCLUS is equivalent to the cone algorithm with R sep = 1. The inclusive dijet cross 
sections vary with the jet definitions. EUCELL corresponds to R sep = 2 and yields the largest 

0 


cross sections. Calculations with the HERWIG Monte Carlo 


reveal that the cross sections 


with PUCELL are larger (almost equal) compared to the KTCLUS cross section for resolved 
(direct) production. The resolved (direct) cross sections are again defined with the x 7 < 0.75 
(x 7 > 0.75) cut. By comparing with the EUCELL ( R sep = 2) cross section we conclude that 
the resolved PUCELL cross section corresponds to a cross section which we would obtain in 
the NLO calculation with R sep ~ 1AR. This is consistent with recent measurements of these 
dijet cross section with different jet definitions |T7|. Here for the enriched resolved 7 sample 
(0.3 < x 1 < 0.75) the cross section for PUCELL is approximately 30% larger than the KTCLUS 
cross section. Both cross sections have appreciable errors, so that definite conclusions about 
equivalent R sep values are difficult to obtain. Furthermore for the measurement of jet shapes 
only inclusive single jets are considered so that it is not justified to read off the exact values for 
R sep from the dijet analysis described above. However, we are confident that for the PUCELL 
algorithm R sep > 1.3 R must be assumed when applied to the calculation of jet shapes where 
the resolved component is dominant. From earlier calculations [3j. Q we know that p depends 
sensitively on R sep . The larger R sep is the broader the jets are. 


The jet shapes for the reaction e + p — > jet + X have been measured for jets with Et > 14 
GeV integrated over four non-overlapping 7 regions (—1 < 77 < 0, 0 < 77 < 1 , 1 < 77 < 1.5, and 
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1.5 < r) < 2) and are presented also in four E r regions (14 GeV < E T < 17 GeV, 17 GeV < 
E t < 21 GeV, 21 GeV < E T < 25 GeV, and E T > 25 GeV) integrated always over the 
same rj interval — 1 < rj < 2. It would not be reasonable to insist that the jet shapes p for 
these different rj and Et regions should be calculated always with the same R sep parameter. 
We expect R sep ~ 1.3 R for the rj interval — 1 < rj < 0 where direct and resolved production 
contribute in nearly equal amounts and it is known that single-jet production in this region is 
dominated by x 7 ~ 1 . Away from the backward (direct) region we expect increasingly broader 
jets the more we approach the forward direction rj —> 2. We know from earlier investigations 
that the single inclusive jet cross section as a function of ?/ for Et > E™ m as measured by ZEUS 
18f does not agree with the NLO prediction in the region 7] > 1 [|T9j, the experimental cross 


sections are larger than the theoretical ones. This excess of the cross section from ij > 1 could 
be simulated in the PYTHIA Monte Carlo |2(J by adding multiple interactions. These multiple 
interactions, which apply only to resolved processes, consist of interactions between partons in 
the proton and photon remnants calculated as LO processes and generated in addition to the 
primary hard scattering. These multiple interactions lead to an energy flow outside the core of 
the jet due to a possible underlying event. In [J7J it was found that these multiple interaction 
effects on the jet shapes are small in the region — 1 < rj < 1 , but increase gradually with rj in the 
region 77 > 1, where an improved description of the data is obtained. In our NLO calculation 
we shall describe this excess of energy outside of the core of the jet with a gradual increase of 
-R se p towards rj — 2 . 


3 Results for Jet Shapes in ep Collisions 


As already introduced in section 1 , for a sample of jets of transverse energy Et, defined with a 
cone radius R, the jet shape p{r, R, Et, rj) is the average fraction of the jet’s transverse energy 
that lies inside an inner cone with radius r < R, which is concentric with the jet defining cone. 
Then the quantity 1 — p stands for the fraction of Et that lies in the cone segment between r 
and R. This is given, up to higher order terms in a s , by 


1 - p(r, R, E t ,j, rj) 


f dE T E T da(e + p —> e' + 3 partons + X)/dE T 
Et,jct(Et,j) lo 


( 8 ) 


where the integral in the numerator is performed over the cone segment between r and R. 
This quantity is calculated from the contributing 2 —> 3 parton subprocesses. For r > 0 the 
integration does not include the collinear singularities which are at r = 0. Therefore it can be 
computed easily from the Et weighted integral of cross sections for e + p —> e ; + 3 partons + X 
over the cone segment between r and R normalized to E t .j times the LO cross section. 1 — p 
in (8) is 0(a s ). 


From ( 8 ) we have calculated p(r, R, Et, rj) as a function of r for the four rapidity intervals 
— 1 < 77 < 0, 0 < r/ < 1, 1 < r] < 1.5, and 1.5 < 7] < 2 taking R = 1 and integrated over 
Et > 14 GeV. Here we identify Rt,j in ( 8 ) with the Et of the jet. For the R sep parameter 
we have made three choices: (i) R sep = 1.4R = 1.4 (since R = 1), which we consider as an 
average value appropriate for the PUCELL algorithm in the ZEUS analysis, (ii) R sep = 2 as the 
maximal possible value, and (iii) an optimized R sep value, which is R sep = 1.3, 1.45, 1.6, and 1.8 
for the four rapidity regions. Our results, plotted in fig. [I], are compared with the preliminary 
data from ZEUS [7[]. We observe, that for fixed R sep the jet profile p does not depend very 
much on 77 except for small r. Furthermore p depends sensitively on R sep . The jets become 
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Figure 1: Jet shape p(r ) for complete single-jet photoproduction integrated over Et >14 GeV 
and four different regions of g. We compare our results using the Snowmass convention with 
R— 1 and three different values of R scp to preliminary 199f data from ZEUS. 



Figure 2: Jet shape p( rj) for complete single-jet photoproduction for r = 0.5 and integrated over 
Et >14 GeV and the same regions of g as in the last figure. We compare our results with one 
variable and two fixed values of R scp to preliminary 1994 data from ZEUS. 
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broader with increasing R sep - Neither R sep = 2 nor R sep = 1.4 agree with the data in the four 
rj regions. However, this is to be expected according to the discussion in the previous section. 
In the first rj interval, —1 < rj < 0, the data are very well accounted for with R sep = 1.3, which 
is the value we anticipated for the ideal situation with no extra broadening due to multiple 
interaction effects with the remnants. With increasing rj the value of R sep must be increased 
monotonically up to R sep = 1.8 in the fourth rj interval, 1.5 < rj < 2. Due to the additional 
broadening from multiple scattering effects towards r/ = 2 we expect this increase of R sep with 
increasing rj. Since we are unable to predict the R sep values for the four rj intervals in advance, 
the optimized R se p ’s come just from the fit. However, we know that 1.3 < R sep < 2.0 in the 
whole rj range. This way, with choosing R sep larger than 1.3, we simulate in addition to the 
modelling of the experimental jet definition broadening effects which are not taken into account 
in the perturbative calculation and which may arise from multiple interactions or other effects. 

If we fix r = 0.5, the dependence of p{r) as a function of rj over the range —0.5 < rj < 1.75 
is shown in fig. ||. For fixed R sep the jet shape p{r = 0.5) depends only marginally on rj in 
contrast to the experimental data [7]. As to be expected, the rj profile is in agreement with the 
data for the optimized R sep (denoted R sep = fit) given by the full curve in fig. || 

In the ZEUS jet shape measurement, p is also presented as a function of r for four non¬ 
overlapping Et regions: 14 GeV < Et < 17 GeV, 17 GeV < Et < 21 GeV, 21 GeV < Et < 
25 GeV, and 25 GeV < Et < 29 GeV, integrated over — 1 < rj < 2. Their results are shown 
in fig. |j compared to theoretical calculations with R sep = 2, R sep = 1.4 and an optimized 
Rsep which varies between 1.4 and 1.5. We see that R sep = 1.4 yields in average a fairly good 
description of the data points over the whole Et region. Only in the first two Et regions we 
need a larger value R sep = 1.5. From this comparison we conclude that R sep increases when Et 
decreases. 

This is seen also when we plot p{r = 0.5) as a function of Et for 15.5 GeV < Et < 27 GeV in 
fig. [| and compare with the data from [^]. Whereas for fixed R sep = 1.4 the jet shape p(r = 0.5) 
decreases only little with decreasing Et the experimental point at Et = 15.5 GeV lies below 
the curve showing that here R sep is larger than 1.4, i.e. the jet broadening with decreasing Et 
IQ is somewhat stronger than predicted for a fixed R sep . 

In fig. §. we show the jet shapes for the four rj intervals and Et > 14 GeV as in fig. |T] with 
optimized R sep for the direct (dashed lines) and resolved processes (dotted lines) separately 
together with the curves, where the direct and resolved contributions are superimposed (full 
lines). It is known that at low Et the resolved component dominates. Therefore the full curves 
lie very near to the resolved curves. The jet shapes for the direct process are, except for the first 
rj region — 1 < rj < 0, narrower than those of the resolved process leading to slightly narrower 
jets when both contributions are superimposed. In the region rj < 0 the role of direct and 
resolved contributions are reversed. Here the direct jet shape is slightly broader than for the 
resolved process. 

The same plots for the jet shapes as a function of r for the four Et intervals as in fig. |3] are 
exhibited in fig. Here all events for — 1 < rj < 2 are included. For p at the smaller r one can 
see that the direct contribution becomes more important with increasing Et as we expect. 






Figure 3: Jet shape p(r) for complete single-jet photoproduction integrated over — 1 < 77 < 2 
and four different regions of E T . We compare our results using the Snowmass convention with 
R = 1 and different values of R sep to preliminary 1994 data from ZEUS. 
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Figure 4: Jet shape p(Et) for complete single-jet photoproduction for r = 0.5 and integrated 
over — 1 < rj < 2 and the same regions of Et as in the last figure. We compare our results with 
one variable and two fixed values of R scp to preliminary 199f data from ZEUS. 
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Figure 5: Jet shape p(r ) for full, direct, and resolved single-jet photoproduction integrated over 
E t >14 GeV and four different regions of rj. We compare our results using the optimal value 
of R S ep preliminary 199f data from ZEUS. 






r r 


Figure 6: Jet shape p(r) for full, direct, and resolved single-jet photoproduction integrated over 
— 1 < 7] < 2 and four different regions of Et- We compare our results using the optimal value 
of R scp to preliminary 199j data from ZEUS. 
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From the comparison in figures |j and |6] it would be difficult to draw any conclusions that the 
direct contribution is needed to explain the data since the general shape of p as a function of r 
is so similar. For this we need more detailed information, for example the q dependence of p for 
fixed r and fixed Et■ This has been calculated in Q for the resolved contribution. The result 
was that p(r = 0.5) for fixed Et and R sep = 2 as a function of rj had a minimum, i.e. a concave 
shape. In fig. j2j. p(r = 0.5), which has contributions from direct and resolved processes, has a 
convex shape with a slight maximum in the considered q range. This occurs for both R sep fixed 
to R S ep = 1-4 and R sep = 2. This different shape of p for fixed r and E T must be characteristic 
for the direct component. This is indeed the case. In fig. we have plotted p(r = 0.5) for 



Figure 7: Jet shape p[rj) for full, direct, and resolved single-jet photoproduction for r = 0.5, 
Et > 14 GeV, and R sep = 1.4. Direct and resolved contributions have different shapes and a 
crossing point near q ~ 0. 

Et > 14 GeV and R sep = 1.4 as a function of q for the direct and the resolved contribution 
separately and also for the superposition of both. As to be expected the resolved contribution 
of p has the concave shape, and the direct contribution is convex at least for rj < 2. The total 
contribution is convex in the region — 1 < q < 2 in agreement with the result in fig. |2j. It would 
be interesting to see the increase of p{r = 0.5) after the minimum at rj ~ 2.0 towards increasing 
rj. In practice this will be difficult, since for increasing rj, R sep must be increased (except at large 
Et), which leads to a decrease of p(r = 0.5) with increasing rj as shown in fig. It is inter¬ 
esting, however, that in the backward direction the shape of p as a function of q is determined 
by the p dependence of the direct component, although the resolved component dominates for 
q > 0. The crossing point is near 77 ~ 0. For q < 0 both components seem to contribute equally. 

From the comparison of the jet shapes as measured by ZEUS with our NLO calculations 
we conclude that the jet shapes can very well be explained with a varying R sep parameter 
describing either additional jet broading effects or/and any possible changes of the merging 
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conditions with respect to the jet merging in the theoretical calculations. In particular we have 
found that for r) < 0, R sep = 1.3 as we expect it for the PUCELL algorithm and _R sep increases 
up to -R sep = 1.8 in the forward direction (see fig. |2|). As a function of E T the parameter R sep 
changes only slightly, when all events for the whole p range are included — 1 < 77 < 2. R sep 
increases in this case with decreasing Et (see fig. [4]) in accordance with the fact that additional 
jet broading effects become more important the lower Et is. In the next section we shall see 
whether these findings can help to describe the jet shapes in pp collisions as measured by the 
CDF and DO collaborations. 


4 Results for Jet Shapes in pp Collisions 

4.1 Comparison with CDF Data 

Quite some time ago the CDF collaboration 0 has presented data for the jet shape p(r, R , Et , rj) 
as a function of r for three E T intervals 40 GeV < E T < 60 GeV, 65 GeV < E T < 90 GeV, 
and 95 GeV < Et < 120 GeV and rapidities in the central region 0.1 < \rj\ < 0.7. The labo¬ 
ratory system at the TEVATRON is the c.m. pp system. These data were compared with the 
NLO results of EKS 0. For E T = 100 GeV, i.e. for the data in the third E t bin they found 
better agreement with the data with R sep = 1.3 instead of the original choice _R sep = 2. The 
measurements by CDF and the NLO calculations were done with R = 1. 

The energy sharing and merging of overlapping clusters and all other jet defining conditions 
in the CDF analysis are the same as described in subsection 2.2 for the ZEUS analysis. This 
means _R sep = 1.3 is the right choice in the theoretical calculation for the treatment of over¬ 
lapping cones. In our calculations we left _R sep to be free to obtain the best description of the 
CDF data. The parton distributions of the proton, which we need as input, are the same as in 
the previous section. Also a s is calculated in the same way as for the ep jets. 

Our results are shown in fig. ^ We present three curves, _R sep = 1.4 and _R sep = 2.0 fixed 
for all Et and a curve with the optimized R sepi which is i? sep = 1.5, 1.4, and 1.3 for the three 
Et bins. Thus _R sep increases with decreasing Et as we expect it from the comparison with 
the ZEUS data. It is clear from fig. that _R sep = 2.0 describes the data in none of the Et 
bins. Second, a small increase of R sep is needed to explain the data in the two lower Et intervals. 

This is observed also when we plot p{r = 0.5) as a function of Et and compare with the 
data J|j] in fig. |[]. Although fixed _R sep = 1.4 is in reasonable agreement, the optimized .R sep 
accounts perfectly for the three data points. We remark that the Et dependence of p at r = 0.5 
in fig. looks very similar to that in fig. |1] for the ZEUS data. Here also the variation of -R sep 
with Et was small (A.R sep = 0.1 similar to AR sep = 0.2 for the CDF data). 

4.2 Comparison with DO Data 

Last year the DO collaboration at the TEVATRON presented their jet shape data || covering a 
similar kinematic range as CDF. Their data were used to populate four non-overlapping jet Et 
ranges of 45-70, 70-105, 105-140, and greater than 140 GeV. The jets were analyzed in a central 
region |?/| < 0.2 and a forward region 2.5 < |?/| < 3.0. In the forward region the data were 
collected only in the first two Et bins. In the experimental analysis jets were reconstructed 
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Figure 8: Jet shape p(r) for single-jet production inpp collisions integrated over 0.1 < \rj\ < 0.7 
and three different regions of E T . We compare our results using the Snowmass convention with 
R = 1 and different values of R sep to CDF data. 
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Figure 9: Jet shape p(Ex) for single-jet production in pp collisions for r = 0.5 and integrated 
over 0.1 < |?/| < 0.7 and the same regions of Et as in the last figure. We compare our results 
with one variable and two fixed values of R sep to CDF data. 
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using a fixed cone algorithm with R — 1. The preclusters were constructed using the Snowmass 
jet direction definitions in (1) and (2). After a stable jet center was found, the jet direction 
was redefined using the DO jet direction definition 

VJ = - ln(tan(0j/2)), (9) 

^ = tan -HY.EyJY.Er,), (10) 

iGjet iGjet 

where the jet polar angle 

Bj = tan-‘( /(JT eR* + (X S K )V X . (11) 

V iGjet zGjet zGjet 

This definition of the jet direction differs from the commonly used Snowmass definition given 
in (1) and (2) above. 

The merging of overlapping jets was done similarly as in the CDF and ZEUS analyses. Two 
jets were merged into one jet if more than 50% of the Et of the jet with the smaller Et was 
contained in the overlap region. If less than 50% of the E t was contained in the overlap region, 
the jets were split into two distinct jets, where the energy in the overlap region was assigned 
to the nearest jet. 


After this procedure the jet directions were recalculated with the DO conventions described 
above. So the main difference in the DO cone algorithm is the definition of the jet direction. 
This influences both the data and the theoretical predictions for the jet shape. It was found, 
however, that the measured jet shapes hardly change, when instead of the DO definition for the 
jet direction the Snowmass definition was used. In the central region r/ < 0.2 the change of p 
was ~ 3% at most (for the inner subcone) and not more than 4% for jets in the forward region 
(2.5 < | )]\ < 3.0) ||. However, on the theoretical side the jet shapes using the two jet definitions 
were much more different. In particular the forward produced jets are much wider with the DO 
definition than with the Snowmass definition. Furthermore the corresponding theoretical jet 
shapes did not agree at all with the experimental data. Thus, whereas the experimental data 
are relatively insensitive to the choice of jet direction definition, the theoretical predictions vary 
appreciably. This latter point was studied further by Glover and Kosower m • They showed 


that DO’s jet direction definition is not as perturbatively stable as the Snowmass definition. 
Therefore they concluded that the DO recombination scheme should not be applied for the 
purpose of making a quantitative comparison with fixed order perturbation theory. Since the 
difference of the measured jet shapes between the two schemes of jet direction definition is so 
small, we shall in the following assume that the data of DO || are obtained with the Snowmass 
definition and shall compare them with our prediction obtained with this algorithm. 


Our results for p as a function of r are plotted in fig. H- separated into the four Et bins and 
for \rj\ < 0.2. The curves in fig. |10| are for fixed R sep = 1.4 and R sep = 2.0 and for an optimized 
R sep (full curve) to produce the best fit to the DO data. The fitted R sep ’s vary between 1.5 
(lowest Et bin) to 1.2 for the Et > 140 GeV measurements. So R sep decreases with increasing 
Et and R sep = 1.3 for the second largest bin agrees quite well with the R sep obtained for the 
CDF data in the same Et range. We observe again, that the variation of R sep is small, Ai? sep 
= 0.3 similar to our result obtained from the comparison with the CDF data. 
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Figure 10: Jet shape p(r ) for single-jet production in pp collisions integrated over \rj\ < 0.2 
arid four different regions of E T . We compare our results using the Snowmass convention with 
R = 1 and different values of R sep to DO data. 




Figure 11: Jet shape p(r) for single-jet production inpp collisions integrated over 2.5 < \g\ < 3.0 
and the two lower regions of E T . We compare our results using the Snowmass convention with 
R = 1 and three different values of R sep to DO data. 
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The same plots for jet production in the forward direction, 2.5 < \r]\ < 3.0, are shown in 
fig- 0- Here we have only two Et intervals. The fitted R sep is R sep = 1.5 in both Et bins, but 
R sep = 1.4 also gives a reasonably good description. 

If we fix r = 0.5, the corresponding p as a function of Et is seen in fig. [12] for the case of 
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Figure 12: Jet shape p(Et) for single-jet production in pp collisions for r = 0.5 and integrated 
over \r)\ < 0.2 and the same regions of Et as in fig. 0. We compare our results with one 
variable and two fixed values of R sep to DO data. 


the central jets. We see that a slight change of R sep is needed to reproduce the data similar 
to the result from fig. |j|. For the DO data this is even more convincing than for the CDF data 
since the range of Et , where experimental data exist, is now larger. 

From this comparison of jet shapes measured in pp collisions at a/s = 1.8 TeV we conclude 
that we obtain a quantitative description of the data when we consider a varying parameter 
Rsep which decreases as a function of Et and is equal to R sep = 1.3 for Et = 100 GeV jets. The 
variation of R se p , needed to explain the measurements in the considered Et range, is rather 
small, AR sep ~ 0.2 in agreement with the findings for jets produced in yp collision at a/s = 
300 GeV. 

So far we have not mentioned that our theoretical predictions depend on the choice of the 
renormalization and factorization scales which we have set equal to Et- Since the jet profile 
1 — p as given in (8) is a LO prediction, i.e. is 0(a s ), the dependence on the scales is stronger 
than in a genuine NLO result which includes corrections O(al). According to (8) the quantity 
1 — p is a ratio of cross sections. Therefore we expect that the dependence on the factorization 
scale cancels to a large extent and 1 — p depends essentially only on the renormalization scale 
p. When p is decreased, a s increases, so that 1 — p increases, i.e. the jets become broader. 
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In principle we could repeat the analysis by varying /i in the interval E T /2 < p < 2 E T as is 
usually done. However, such a variation of p has very little influence on p{r) for r > 0.5 SI- 
For r < 0.5 the change of p(r) is noticeable. Thus the fits to the data might even improve 
in this range of r, when R sep and p are optimized simultaneously. This would not change the 
trend of the R sep variation with decreasing Et- Only the optimized R sep values would change 
slightly. Furthermore, for small r < 0.5, the jet shape p(r) is not predicted very well, since in 
the region r —> 0 higher order corrections become more and more important. The change of 
p with R sep occurs also for large r, where the predictions are more reliable and the data have 
smaller errors. From this we conclude that the R sep values deduced from the data should not 
be very sensitive to the choice of the renormalization scale. The same conclusion can be drawn 
concerning our choice of the one-loop formula for a s versus the two-loop formula fixing the A 
value through the parton distribution of the proton. Making other choices is just equivalent to 
small changes of the renormalization scale. 


5 Summary 

We have calculated jet shapes in photoproduction superimposing direct and resolved contri¬ 
butions and in pp collisions in next-to-leading order of perturbation theory. The dependences 
of the jet shapes on transverse energy, rapidity, and inner cone extension were compared with 
recent data from ZEUS at HERA and from CDF and DO at the TEVATRON. The numerical 
results were obtained using the Snowmass cone algorithm yielding good agreement with the 
data, if the problem of merging overlapping jets and additional jet broadening effects are taken 
into account by varying the parameter R sep as a function of transverse energy and rapidity. 
The optimal values and variations of R sep with Et were found to be rather similar in low Q 2 ep 
and in pp collisions. Furthermore, it was shown that the convex shape of the jet profile rapidity 
distribution in photoproduction is due to the direct photon contribution. 
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